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This report is an attempt to investigate the influence of dipole-dipolc coupling in a 
paramagnetic spin system at low temperatures. It consists of two parts. The first i)art is 
a discussion of the use of Cm=Ch^o for a system with mutual interaction. It is pointed 
out that only if the external field is large compared to the internal field is this equation correct. 

The other part consists of a calculation of higher order correction of the di])ole-dipole 
interaction on a system of paramagnetic spins which is subject to a crystalline field which 
w^e chose of the y2,o type. The total Hamiltonian consists of a spin Hamiltonian in accord- 
ance with this symmetry, a term representing the external magnetic field and the dipole- 
dipole interaction between the spins. The partition function is calculated by means of the 
Schwinger trace formula considering a representation in which the first two terms of the 
Hamiltonian are diagonal. The trace of the density matrix can be expressed as the trace of 
a product, one factor is the density matrix of the noninteracting spins, the otlier factor consists 
of a sum of comnuitators. Tlu^se commutators are workc^d out in detail and the result is 
given in the form of a finite series over the quantum number m. There seems to be; no ob- 
vious way to perform these summations. 



1 . Introduction 

In a system of paramagnetic dipoles in a crystalline field the magnetic specific heat at con- 
stant magnetization is equal to the specific heat in zero magnetic field. This is easily seen from 

and the fact that the magnetization in a system without any spatial correlation (i.e., ideal gas- 
like) is always proportional to H for small fields, hence M as well as {dMjdT) will go to zero if 
there is no field and we have Cm=Ch=^Q' 

As a consequence of the lack of correlation M=0 if H=0, hence the relation should be 
Cm=q=Ch=Q' 'I'he general question when Cm=Ch=q is discussed in appendix 1. 

If there is no crystalline field both these quantities are zero and the magnetic specific heat 
at constant Tly^O is merely determined by the Curie law 

C«=^ (1.2) 

where C is the Curie constant. 

For the case of dipole-dipole interaction between the spins the value of the specific Iieat at 
constant magnetic field was calculated by Van Vleck in 1937 [1].^ The result can be expressed 
in the limit of high temperatures in the form of an internal magnetic field Hi, which has an average 
value zero, but a root mean square value H^, that is non-zero; then 



Figures in brackets indicate the literature references on page 118. 
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We want to calculate Cm again, the reason being that the ratio of Cm to Ch determines 
among other things the real part of the susceptibility at high frequencies. The result that 
Cm=Gh=q was used in the calculation of this ratio [2], although this is correct only if a Curie law 
holds. Although the error is not very important, at least not at external fields large compared to 
the internal field, there is a certain inconsistency in the fact that Hf^^O is in contradiction 
with the Curie law, which assumes the total absence of any mutual interaction. 

In order to reconsider the equation of state, one can draw up a modification of the Curie 
equation by replacing the field H by either H-\-Hi or by ^JH^-^Hi and then recalculating the 
difference between Ch and Cm again; this is not equal to zero ior H=0. As a matter of fact it 
is infinitely large with the second choice 

jif-cg±g.'->c„-c:„— ^<g+g->' (1.4) 

The failure of such attempts is due to the fact that the Curie law is a high-temperature 
expansion in the case where the internal field is small. That is, the internal field itself is 
temperature dependent in such a way that it will vanish for T-^oo . The difference between Cm 
and Cff is now 

since for high temperature Hj is small and we find that this difference goes to zero for H=0 and 
that the correction is proportional to T~^. This could have been obtained in a much simpler 
way if one assumed the Curie-Weiss law for the magnetization 

M=Cj^' (1.6) 

which leads to 

T 

Ch—Cm^ — CH ^rji_0Y' (1-7) 

The introduction of interaction between the spins is not in contradiction with the assumption 
that AC goes to zero for H—>0. 

The considerations above were mainly intended to provide a reasonable excuse to introduce 
an internal field by means of those terms in (kT)~'^ which arise from the dipole-dipole interaction. 
These terms can be classified in two groups, classical and quantum mechanical. The quantum 
mechanical terms are the result of the noncommuting of the operators. For zero external 
field they occur only if there is a crystalline field splitting. It is the main purpose of this work to 
calculate these quantum mechanical effects in the internal field. The method followed is straight- 
forward. The effects should be observable in the '^intermediate" temperature region, that is, 
the high temperature expansion should hold far enough down, but not too close to the Curie 
temperature of the dipole-dipole interaction [3]. 

2. Operator Expansion 

We will separate the Hamiltonian into two parts: the first part, ^i=a/kT, contains the 
diagonal terms, the spin Hamiltonian in the crystalline and the external field; the second part, 
g^2=b/kT is the dipole-dipole interaction. It has no diagonal terms. 

The expansion of the exponential needed for the evaluation of the partition function is 

g-(a+&)^g-a_^^ (_l)n^~a^^^^ (2.1) 

n = l 
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As is well-known [compare 4] Sn contains the 6-operator r^ times. If one actually wants to 
have concrete values for Sn one lias to perform another expansion, this time with respect to 
a. The practical question is how far do we calculate this double power series. The term 
e"°.6 contained in Si is a product of a purely diagonal term and a purely nondiagonal term. 
Hence the trace is zero [5], which means that the average value of the field created by the 
dipole-dipole interaction at a certain point is zero. The field has, however, a mean square 
value, and the first term in ^2 is the major contribution to its value. This is a purely '^classicaP' 
term in the sense that it is also nonzero for conmmting operators. All first terms in S^ are 
classical, but since we assume that the dipole-dipole interaction is small compared to the 
thermal energy of the system the description on the basis of the first term in S2 will be sufficient 
for the classical part as long as there is no evidence of an approaching cooperative phenomenon. 
Another argument for dropping the next terms is to notice that the range of the dipole-dipole 
interaction squared is r~^, and all higher classical terms with even 7?.-value will give rise to weak 
interactions of extremely short range unless the temperature is so low that a strong correlation 
between the spins occurs. 

The expressions for Si and S2 are given by 

Si= f r/s,6(.s,); S2- f <l.s, r\/s,b(.s,)h(.s.^ (2.2) 

Jo Jo Jo 

and expanding b (n) : 

b(s)=e"'be-''' = b + 4a, b] + . . . (2.3) 

The result of the integration is 

+ m(i, b] + ^^b\a\a, b\\ + l\a, b]b + l\a, bV+l[a[a, b]\b. (2.4) 

None of the terms of S^i contribute to the trace since they are the prochict of a number of 
diagonal matrices and one nondiagonal matrix. 

The terms of S2 up to the second oi-cUm' in b are 

lb'+lh[a,b]+l[a,h]b. (2.5) 

It is shown in the app. 2 tliat the only contribution to the weighted trace is 

tr {e-Uib'+\\b\a,b\\)]. (2.6) 

3. Calculation of the Double Commutator 
In order to calculate the double commutator we introduce the following notation 

^2= ciV 12 '^^'S''-3(^'^^")(^'^>')/^» = ^ y-, S'aQ'4Sl. (3. 1) 

k?^l >"lcl a, /3 kl£l 

Where the position indices are written as superscripts and the matrix 6 is given by: 

0ii=!7'/S'(2/'+s'-2x^)/r^; Q^^^Q^.^ifff'xy 

Q^z=iffi\x'+f-2z')lr'; Qi,=eu = (f^'xz (3.1a) 

leaving out all superscripts. 

Since the leading part of the Hamiltonian is of the form (a, b, and c are constants) 

^1= S aS', + bS,+c (3.2) 

positions 
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we have to calculate the commutators C„ and double commutators Z>„: 



(3.3) 



where all spins are at the same location. In order to avoid confusion with the superscripts^ 
which have been omitted, the powers are placed outside parentheses. Table 1 shows that 
there are four nonzero values of Di. However two have zeroes in the diagonal, hence if multi- 
plied by the (diagonal) weight factor only a=P=x or y will give a contribution. The same 
holds for D2. The trace of the two elements of C is fi^m. The trace of the two elements of 
D is —fi^[3m^—S(S-\-l)]. This has to be inserted in the multiple sum over all positions 
containing two Q's and multipled by the weight factor exp {SKiJkT). 

Since both Hamiltonians are series over all positions and since all commutators with 
respect to two different positions are zero we have to figure out in which ways these two 
summations can ''overlap'\ The first commutator gives (k^^l): 

[^2, ^i]=«x: fs si@^',s'„ (-s?)1+6S rS'S'^e"^'^^, is?y']- (3.4) 

klm ^ a^ -J klm L «^ ■ ^ J 

The summation over m contains only two terms m=k and m = l. These terms are equal 
since the matrix is symmetric in both the upper and lower indices and since ^ 5^ Z the operators 
commute. 



k9^l a,0 !C9^l off 



(3.5) 



The double commutator introduces a sum over four positions, say k, l{k^l) and i, j{i^j). 
There are six different terms (table 2). The first four give all the same contribution because 
of the symmetry of mentioned above, and so do the last two. Writing for short 

we have 

oifi \ jkl 
a'fi' all different 

+lS'a', [S'a, f{S')]mgQ'J>^'S'^,S',-[Sl,, [SI, f{8')]WJgQ»,^,[S'^,,S',] } ), (3.6) 



Table 1 





Di/m' 
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-Sz 
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-S^ 
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Sx 
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Diliifi)^ 



fi=x 



(S^Sy+SyS^) 



-{SySxi-SzSy) 

2{S!-Sl) 

(S^Sy+SyS^) 



Table 2 
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-,6 
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116 



iising- for the second set of tenns the following commutator identity: 

[AB, CD]=AC[BD]-\-[AC]BD-[AC][BD] 
if only 

[AC]y^{) i\m\ [BD]^0 

and all other conunutators are zei'o. Only the diagonal part of this double commutator 
contributes to the partition smn, and as a ]'esult of this tlio suinmation over «, /3, ex and ^' 
has a number of restrictions: 

First term: 0L = a =x or y; P=l3^ = z. 

Second term is always zero since [S^, f^So)] has no diagonal elements for every p\ 

Third term: a, a =x or y, and /3, fi'=x, y or ?/, x. 

Fourth term: a = a=x or y and /3', fi = z. In case it is x,y it will be canceled by the 
diagonal contribution from y, x. 

Fifth term: a = a =x or y and /3', /3=x, y or y^x. 
The diagonal parts of the spin operators to the second power are: 

fiS 
and using table 1 we firid (condoning the first and fourth term): 

ky^l 

+hS'immi+Q'yi'^M)siS'\+T: 4 2: ([si,j{S'^)\si,(dim''Ai^',-,s'A (3.7) 

The last term is calculated with the lielp of 

[^V^, Sl\S^ = {T2llS,^ri')S^S^ + hKS^-?>S'l)±h'S,T'2hS,^^^^ (3.8) 

The summation over a, 13, a , 13^ contains 8 terms. It turns out, however, that 4 cancel mutu- 
ally: if a and /3 are equal (both x or both ?/) the subscripts a^' are either x and y or y and x. 
Both have the same coefficient QxxQx^y l^ut tlie comnmtators have the opposite sign. The 
fact that these cancel could have been seen on the following general grounds. A comnmtator 
is anti-Hermitian and a double commutator is Hermitian. Hence the last factor has imaginary 
diagonal elements and the total expression has real diagonal clenients, and the first factor must 
have an imaginary diagonal part: 

is„,f{s,y]s,=-[s,j(s,)]s, 

=l({S+J(S.)]S--[S-,JiS,)]S+ 

= -^ [a{2hS,iS'-Sl)-ri'S,} +bh{S'-Sl)]. (3.9) 

The final result for tlie last term is: 

4 2 h'l(.e'.'M'y)-(QVyy] [2aSKS'-S!-hh^)+bS,{S'-Sm. (3.10) 

k9^l 

In order to obtain a value for the partition functions the trace has to be calculated nu- 
merically since the sums of the type 

J2p= S mV^'+'^^^ (2>=0,1,2,3) (3.11) 

m=—S 

cannot be evaluated in closed form. 

117 



4. Conclusion 

The complication of the calculation lies mainly in the fact that the commutator of two 
spin operators is not a c-number but again an operator. This makes many tricks used to 
evaluate exponential operators useless. For instance, one could have written the original 
total Hamiltonian with step up and step down operators S'^ and S~. The expansion of the 
exponential would have created a polynomial containing S+, S" and Sz to certain powers. 
A general conclusion is that only those terms which contain an equal number of S'^ and S~ 
at the same location will have a nonzero diagonal element. To calculate the actual value, 
however, one has to shift the S'^ and S~ operator such that they are all at the end or the be- 
ginning of the monomial and in conjunction to each other. This rearrangement of operations 
is accomplished with the help of Wick^s theorem if we deal with creation and annihilation 
operations. It is still possible to formulate a Wick-like theorem for spin one-half. This 
was done by Kenan [6]. It introduces an additional numerical factor into each diagram. 
For spin not equal to one-half, his treatment fails. The crucial point is that only for spin 
one-half is (Sz)^ a c-number. Otherwise Sz is an operator and hence can not be taken outside 
the trace and the contraction operator cannot be defined. These considerations indicate 
that it seems unavoidable to calculate the trace straightforwardly. Actually the straight- 
forward calculation consists of doing the same steps as in the Wick- theorem. 

Finally, the summation of the positions has to be performed. It is possible to restrict this 
to the immediate environment since the radial dependence decreases very rapidly. 
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6. Appendix 1 

We want to investigate under what conditions Cm=Oh=o for arbitrary H. To calculate 
the specific heat at constant H we find from the Maxwell relation dS/c>H=dM/(>T 



the second term at the right-hand side of (11) is 
and the condition Cm^Ch^o gives 

ns 



wliicli leads to the differential ec|ualion 

KdHj dT' dTdH dT c>H^\dTj dW 
which is fulfilled by 

where / is an arbitrary function and Ci and c-z two constants. The initial concHtion forces 
us to take Ci^O. For kT larger than the crystal field splitting this condition is fulfilled. 

7. Appendix 2 

From the sum and difference 

b[a, b] + [a, b]b=ab'-¥a b[a, b]-[a, b]b^\b, \a, b\\ 
we find 

h[(h h\ = h[b\a, b\\ + \a¥-Wa [a, b]b = h[b, [a, b\\-lab'+Wa. 

Hence the relevant terms in (2.5) are 

lb[a, b\ + \[a, b]b=-^,\b, [a, b\]+^\jab'-^,b'a+i[b, [a, b]\-iab''+Wa=l\b, \a, b]]-^ab'+^^b'a. 

Since a is diagonal we have 

ii:[U(a)aP)} =tr. iUiaWa)}, 

lience the contribution to the trace is 

iv.{Mm<',b]]}. 

(Paper 6SA1-259) 
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